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We present an extensive study on inhomogeneous chiral condensates in QCD at finite density in 
the chiral limit using a generalized Ginzburg-Landau (GL) approach. Performing analyses on higher 
harmonics of one-dimensionally (ID) modulated condensates, we numerically confirm the previous 
claim that the solitonic chiral condensate characterized by Jacobi's elliptic function is the most 
favorable structure in ID modulations. We then investigate the possibility of realization of several 
multidimensional modulations within the same framework. We also study the phase structure far 
away from the tricritical point by extending the GL functional expanded up to the eighth order 
in the order parameter and its spatial derivative. On the same basis, we explore a new regime in 
the extended GL parameter space and find that the Lifshitz point is the point where five critical 
lines meet at once. In particular, the existence of an intriguing triple point is demonstrated, and its 
trajectory consists of one of those critical lines. 
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I. INTRODUCTION 

How matter behaves as a whole under extreme con- 
ditions, such as high temperature and/or density is one 
of the central questions in hadron physics that, in prin- 
ciple, should be answered by quantum chromodynamics 
(QCD). At high temperature QCD unambiguously pre- 
dicts a new, extreme form of matter called the quark- 
gluon plasma, and its properties are now under exten- 
sive investigation from both theoretical and experimental 
sides, the latter facilitated by heavy ion colliding facili- 
ties such as RHIC and the LHC. On the theoretical side, 
the most powerful tool to investigate this region is the 
lattice QCD simulations which makes possible ah initio 
calculation of thermodynamic quantities at high temper- 
ature. 

On the other hand, at finite density lattice simulations 
still remain less predictive because of a serious problem 
due to the complex fermion determinant. Our knowledge 
of the forms of matter at finite density is limited except 
for ultimately a high density regime where we know, from 
pcrturbativc QCD computations, matter forms a color 
superconductor of the color-flavor locked type. Even 
though the form of matter at accessible density is still 
obscure, we might expect to have two kinds of phase tran- 
sition (or crossover) when ordinary nuclei is compressed 
by an external gravitational pressure: One is the dccon- 
finement [l| and the other is the chiral restoration Q 
which is the main topic of this paper. 
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QCD vacuum is considered as a sea of the chiral con- 
densate 0] , namely, qq condensate which may arise due to 
a strong correlation between a quark and antiquark. The 
assumption has proven to be useful for interpretation of 
hadron spectroscopy as well as hadron scatterings at low 
energy, yet it is only recently that its formation in vac- 
uum was directly checked by lattice simulations of QCD 
3. When the baryon chemical potential is introduced, 
it breaks the charge conjugation symmetry, bringing an 
imbalance, i.e., a net excess of quarks over antiquarks. 
Thus it is possible for the chiral condensate to be broken 
by such imbalance. In the model calculation this actually 
takes place through the suppression of qq formation by 
the Pauli blocking 

How the entire chiral restoration takes place is, how- 
ever, still under the veil. Most of the Nambu-Jona- 
Lasinio (NJL)-type models predict the first-order phase 
transition. But it is possible that this transition proceeds 
hierarchically via several intermediate steps, such as me- 
son condensates. In particular, possibilities of inhomo- 
geneous chiral condensate [chiral density waves (CDW)] 
have been attracting a lot of interest recently |a-[lO|. 
So far proposed CDWs can be classified into two major 
groups: one is the Fludc-Fcrrell (FF)-tvpe \VV{ and the 
other is the Larkin-Ovchinnikov (LO) type [i2|. CDWs in 
the former group are characterized by some anisotropic 
pairing in momentum space, while those in the latter 
group are described by a characteristic modulation in 
real space. 

CDWs can be regarded as fcmtoscalc phase separa- 
tions as a consequence of the competition between the 
chiral condensate a = {qq) and quark chemical potential 
/i; the latter drives a stress that promotes the population 
of quarks, n = (q^q), that is, an imbalance in the sea of 
quark-antiquark condensate. In the case of the FF-type 
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CDWs, quarks as impurities are accumulated in some 
favorable directions on the Fermi surface where the pair- 
ing is absent, while in the LO-type CDWs the condensate 
modulates in real space, taking nodes where quarks are 
populated without any friction [T3| . In a nutshell, it can 
be said that the CDW is a halfway state produced by a 
compromise between two conflicting effects: the qq con- 
densate and the quark occupation q^q. 

Essentially the same picture applies to conventional 
FFLO superconductors in condensed matter physics [ll|, 
[12} . It was recently reported that a kind of such exotic su- 
perconductors was observed experimentally in 2D films 
of a heavy fermion compound, CeCoIn5 [IJ]. In these 
systems, an external magnetic field h drives the Pauli's 
paramagnetism via yielding a mismatch in the Fermi mo- 
menta of two species, in this case, spin up and down. This 
paramagnetic stress competes with the spin singlet pair- 
ing, and as a consequence, the FFLO state may appear as 
a halfway state between a complete paramagnetic state 
and a uniform pairing state; when h exceeds some critical 
value, the spin density n3(x) = {pr^f) = ?T-t(x) — 
starts to accumulate near the points where spin singlet 
pair density A(x) = {f^{iT2)f) = (ti — it) takes nodes 
in the real (or momentum) space. 

Possible formation of inhomogeneous condensates has 
been discussed in various contexts of hi gh e nergy physics. 
These include neutral pion condensates |l5l|. charg ed pion 
condensates [l^, color superconductors [l7|42 ll chiral 
magnetic spiral in a strong magnetic field [22|, crystal 
phases in (1 -I- l)-dimensional QCD in the large limit 
p3l |. and spiral phases in the quarkyonic phase [23 - l26| . 
Inhomogeneous chiral condensates have also been stud- 
ied extensively in the Gross-Neveau (GN) model with 
a discrete chiral symmetry [27l - [3l| and also in the cor- 
responding model with a continuous chiral symmetry, 
i.e. the (1 -I- l)-dimensional NJL model in the large N 
limit j32| . 

Even if the crystal structure is restricted to ID modu- 
lations, its impact in the QCD phase diagram is drastic as 
demonstrated in [l^] ; it was shown that there is a domain 
in the phase diagram where the chiral condensate devel- 
ops a crystallization in a solitonic form characterized by 
the elliptic function. Accordingly, the tricritical point 
turns into the Lifshitz point where three phases, a sym- 
metric (Wigner) phase, an inhomogeneous phase, and 
a homogeneous symmetry-broken phase, meet at once. 
The analyses were based on the Ginzburg-Landau (GL) 
functional expanded in the chiral order parameter and its 
spatial derivatives up to sixth order, which offers a mini- 
mal model-independent description of the QCD tricritical 
point in the chiral limit (33| . Despite the development in 
understanding the ID modulated chiral condensate, only 
a few works have been devoted to exploring multidimen- 
sional modulations [2^, [s^l • 

The aim of the present paper is twofold: first to seek 
the most favorable crystallization pattern and, second, 
to explore the phase structure away from the tricritical 
point. In connection with the first point, we note that 



even in the GL framework it is not so evident that the 
solitonic chiral crystal characterized by Jacobi's elliptic 
function gives the absolute ground state. What is shown 
is that it certainly constitutes a sufficient, particular so- 
lution at sixth order of the GL functional. We thus exam- 
ine other crystal structures as candidates of the ground 
state, including the possibility of multidimensional crys- 
tal structures. For a ID structure, we search for the most 
favorable structure taking the most general Ansatz, the 
condensate expanded in the harmonic series [2l|, [s^ . As 
for a higher dimensional structure, we only examine sev- 
eral specific Ansdtze. 

In order to accomplish the second purpose, wc ex- 
tend the existing GL functional up to eighth order in 
the condensate and its spatial derivative. This is be- 
cause the GL functional at sixth order provides us with 
a minimal description of the critical point, where it fol- 
lows from a simple scaling rule that the phase structure 
stays the same even if the region close to the critical 
point is zoomed out. Moreover there is one more pos- 
itive reasoning for this extension: It was shown in the 
context of a superconductor under an external magnetic 
field that such higher order terms should be taken care 
to correctly describe phases off the tricritical point (36l |. 
Starting from the BCS model and making use of a quasi- 
classical (Eilenbcrgcr) approximation to the Bogoliubov- 
de Gcnne equation, they numerically found that another 
type of FFLO phase named FFLO-II shows up in the 
phase diagram; it replaces a large part of the inhomo- 
geneous phase, which brings a new critical end point lo- 
cated off the tricritical point. In order to understand 
their results within the GL framework, one needs to work 
at least at eighth order as pointed out in the paper. 

The paper is organized as follows. Section |TT] is mostly 
devoted to a review on the GL treatment of inhomoge- 
neous chiral condensates made in [2l| . We make physics 
behind a formation of inhomogeneous condensate trans- 
parent. In Sec. mil we examine several Ansdtze for mul- 
tidimensional crystal structures, as well as the most gen- 
eral Ansatz for ID modulation. We finally address in 
Sec. llVl the question of how the phase structure expands 
at large distance from the Lifshitz point by extending the 
previous GL functional up to eighth order. We summa- 
rize the present paper in Sec. |Vl 



II. GENERALIZED GINZBURG-LANDAU 
APPROACH ON INHOMOGENEOUS PHASES 



Nickel studied the possibility of inhomogeneous chiral 
phases in QCD using a generalized GL approach [loj . 
Here we review the framework. In order to describe the 
phase structure close to the tricritical point, we need to 
retain up to sixth order in the chiral order parameter in 
the GL Lagrangian. Following [l^l , we extend the GL La- 
grangian such that the theory allows an inhomogeneous 
phase. 
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A. Generalized Ginzburg-Landau functional 

To incorporate the energy gain/loss with respect to 
space variations of the order parameter, we employ the 
gradient expansion. A generahzed GL functional in the 
chiral limit up to sixth order in the chiral order parameter 
and its derivative is given by 



original GL parameter space {02, 0:4, a^} according to 
the scaling relation between a2 and 772 in Eq. ([3]). 

In the following we use {0:2, M, x} instead of {a2o:e, 
My^ctQ, x/y/og} to avoid notational complications. They 
should be understood as dimensionless ones measured in 



the proper units, i.e., [ag ], [«( 



and 



l/2i 



;(A/(x)) = ^M{x)^ + ^M{x) 



aib 
4 



(VAf(x))2 



(1) 

Nickel derived, starting from an NJL model, relations 
among parameters {«„} up to total derivatives: Q!4b — 
04, aeb = Sag and aec = ae/2. Hereafter, we treat these 
three parameters {a2,a4,ae} as independent GL cou- 
plings. In the NJL model, only quark loops contribute 
to the GL couplings, which may provide a good descrip- 
tion at high density regime. Then the generalized GL 
functional to this order becomes 

w(Af(x)) = ^M{^y + ^{M{x)^ + {VMf) 
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(A/(x)6 + 5M2(VM)2 + i(AM)2) 



(2) 

We can then analyze the phase structure of quark mat- 
ter in the GL parameter space spanned by {a2, a4, a^}, 
which can in principle be mapped onto the phase diagram 
in the (/i, T) space in real QCD. 



C. Condition for soliton/condensate formation 

For a4 > the phase structure is simple. One has 
a uniform chiral symmetry-broken phase (M ^ 0, x^B 
phase) for 7^2 < and symmetry restored phase (M = 0, 
Wigner phase) for 772 > 0, separated by a second-order 
phase transition. On the other hand, phases for < 
are rather rich [l^]. When the condensate is restricted 
to be constant in space, there are only two phases: the 
xSB phase at small 7/2 and the Wigner phase at large 
772, separated by a first-order phase transition at 772 = 
3/16 = 7^2. Once the chiral condensate is allowed to 
have variations in space, the first-order phase transition 
splits into two second-order phase transitions, one at 772 = 
5/36 < 772, the other at 772^ = 3/8 > 779, and in between 
an inhomogeneous phase shows up [lO| as we will discuss 
below. 

Let us start with a one-dimensional sinusoidal chiral 
density wave @, 0] (the LO state) characterized by 



AfLo(x) = \/2A/iisin(fciiz). 



(5) 



B. Dimensional and scaling analysis 

Parameter ae should always be positive for the ther- 
modynamic stability and it has dimension A~^ with A 
being some energy scale. Then we use ae to specify the 
unit of energy scale. Moreover we can use {a^l to ad- 
just the magnitudes of the fourth- and sixth-order co- 
efficients. To be precise, let us introduce dimensionless 
variables {772, w, 7n, x, V} as follows: 



a2 



X 



V2 



M ■ 



„l/2 



(3) 
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Then, in these dimensionless bases, the GL functional 
takes the form 



(4) 



Now we see that the problem has only one continuous 
parameter 772. What we have to work out now is to find 
the phase structure in 772 space both for positive and 
negative cases, and then map them back onto the 



Averaging the GL free energy density over the Wigner- 
Seitz cell, optimizing the wave vector kn with respect to 
Mil, and finally expanding the GL free energy in Mn, we 
obtain for < 0, up to quartic order in Mn, 



O = (cj(A/lo(x)))ws 



(f 



16; 



4^'^II 



0«) 



(6) 



From this, it is clear that the system undergoes a second- 
order transition from the Wigner phase to the LO state 
when the parameter 772 crosses 772^ = 3/8 from above. 

Now we move on to the region of the lower value of 7/2, 
where the chiral symmetry is broken by a homogeneous 
condensate. Consider the situation where a single soli- 
ton (a kink in ID [23, HO, [Ml , actually a domain wall in 
3D [lOl) is formed in a sea of homogeneous condensate. 
When 772 is small, there is an energy cost for the creation 
of such a kink in the homogeneous condensate. However, 
this energy cost decreases with increasing 772 , and at some 
critical point which we denote by ri\ it vanishes. To see 
this, we first set an Ansatz for the domain wall profile as 



A/ss(^) = A/itanh(A:iz), 



(7) 



where Mi is the homogeneous condensate, being a func- 
tion of 7^2. The energy cost accumulated within the do- 
main wall surface per unit area can be defined by 

me{Mi,h;Ti2) = /5//2dz(w(A//ss(z))-w(Afi)). (8) 
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where £ is the system length in the z direction which 
will be taken to be infinity. Optimizing the value 
of ki requiring that the energy cost F{ri2-,ki) — 
\mii^oo[iSfli{Mi,ki;ri2)] takes minimum, we obtain the 
relation h = Mj. The function f{7]2) = ^(772,^^1(772)) 
is depleted in Fig. [TJa), from which we see that it is a 
monotonieally decreasing function of 772, and at the crit- 
ical point which we denote by 772 it crosses zero. At this 
point, a soliton will be formed spontaneously without 
any energy cost. It is easy to derive 772 = 5/36 where 
Mi = y^5/6, and in fact 7(772) is expanded in powers of 
A = (772 - 5/36), as /(772) = -|A + fA^ + 0{A^). In 
the figure, we denote by a sohd circle 772 = 3/16, and 
by the cross 772' = 2?72 for comparison. Obviously the 
domain wall formation is faster to take place than the 
first-order phase transition to the Wigner phase when 772 
is increased. We conclude that some kind of inhomoge- 
neous phase should develop in the window 772 < 772 < 772'. 

In order to elucidate on what drives the formation of 
a domain wall in the spatially uniform x^B phase at the 
critical point 7^2, we show in Fig. [Ijb) the energy den- 
sity profile of a single soliton by the bold line. The spa- 
tial profile of the soliton is also shown by the thin line. 
Several contributions to the energy density are shown 
separately: The dotted (red) line shows the energy loss 
through the homogeneous parts in the GL functional, 
while the dashed (blue) and dot-dashed (green) lines 
show the contribution from the derivative terms at the 
fourth and sixth orders, respectively. We confirm that 
the derivative term at the fourth order is responsible for 
the formation of a kink in the chiral condensate. 



D. Solitonic solution at sixth order of GL 
expansion 

We have seen that, for 014 < 0, the first-order chiral 
restoration at a2 = 7/2 0^4 /ag is replaced by two second- 
order critical lines, one at a2 = 772a|/a6 representing 
the onset of the formation of a single soliton in a ho- 
mogeneous sea of chiral condensate, and the other at 
ck2 = ii¥'^4:/'^6j corresponding to the chiral restoration 
from an inhomogeneous chiral phase. In short, the chi- 
ral restoration becomes a smooth, hierarchical two-step 
transition. So now the question is how does the system 
develop from a single soliton state to a sinusoidal LO-type 
state as 772 is increased from 77^ to 772^? This was already 
addressed in the paper [3l|, [STllSSj where it is shown that 
the solution which covers these two situations is the soli- 
tonic condensate characterized by Jacobi's elliptic func- 
tion. Here we closely follow the discussion in [sj and 
demonstrate that, if the modulation is restricted to ID, 
Jacobi's elliptic function gives in fact an adequate solu- 
tion to the problem. Let us start with the Euler-Lagrange 
(EL) equation for the inhomogeneous chiral condensate, 

[ dyw(M(y)) = 0. 



(a) 

0.15 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
12 



(b) 




z 



FIG. 1. (color online), (a) Energy per unit area required for 
the creation of a single soliton as a function of 772. The solid 
circle and the cross placed on the 772 axis represent points 
7?2 ~ rjl and 772 = 772. See text for details, (b) A single 
soliton profile Mss = Mj tanh(Mi2;) with Mi = a/5/6 being 
the mass at 772 — 1^2 a-nd the associated energy density pro- 
file. Contributions from homogeneous parts, and those from 
derivative terms (at the fourth and sixth orders separately), 
are also shown. 



The condition results in the following fourth-order ordi- 
nary, but nonlinear differential equation (when Q!4 < 0): 

= A/(4)(2) + 3M"-10[M(Af')^ + A^^^^"] 
+6772M - 6M3 -\- 6M^. 

We need to solve this equation to find a suitable ID mod- 
ulation. There is no systematic way to find out the gen- 
eral solution to a fourth-order nonlinear differential equa- 
tion. However it is not so difficult to see that Jacobi's 
elliptic function is a particular solution to Eq. © [37| . 
Let sn(z, v) be Jacobi's elliptic function with v being the 
elliptic modulus. We then set 

Msn{z)^ Moi^sn{kz,iy), (10) 
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where k and Mq are constants0. Following the discussion 
in [37l |. we can show that this function obeys a fourth- 
order differential equation 

-^(12 - B) [MsniM^nf + M^M^'y 

+ '^^){3-B + A) + 1| (4 - B) M^^, 

(11) 

with j4 and B being arbitrary numbers. Comparing 
Eq. (|TO|l and Eq. ((TT|) . we see that A/gn constitutes a 
sufficient solution to the EL equation (|10p when the fol- 
lowing five algebraic equations are all satisfied: 

6?72 = [A{v^ + 1)2 - Bv^] , (12a) 
-3= -(l + A)fc2(,.2 + i), (12b) 

10=^(12-5), (12c) 
-6^-^{l + iy^){'S-B + A), (12d) 

6=||(4-i?). (12e) 

From Eqs. (fT2c)l and (fT2e)) . we see {Mo/k,B) = (1,2) or 
{Mo/k,B) = (2,-28), but the latter results in complex 
fc^ and J^. Then we take the choice {Mo/k,B) = (1,2), 
and in this case we see Eqs. (|12bp and (jl2dp are degener- 
ate so that we are left with only two algebraic equations 
p2ap and (|12bp . whereas we still have three unknown 
parameters, v, fc, and A. From these two equations we 
can solve v and /c as a function of the continuous param- 
eter A (and 772). Denoting these functions as kA and i'a, 
we arrive at a one-parameter solution group to the EL 
equation: 

Msniz;A) = kAVAsn{kAZ,iyA)- (13) 

We stress that for any value of A, as long as < z^a < 1 
and kA> are both satisfied, Msn{z; A) gives a solution 
to the EL equation. The potential n{Msn{z] A)) depends 
on A so that we need to look for an optimized value of 
A demanding that it takes a minimum. This is easy to 
work out numerically for rj\ <rj2 < rf^ ■ We see A — > 4/5 
as 1]2 ^ +O1 where va approaches unity [va — > 1) and 
kA ^ ki = -^/S/G, 

lim Msn(z;^) = A:itanh(fcia;), (14) 

while A — > 1 as the other end 772 = rj^ — Q is approached, 
where va goes vanishing and kA — > fcii ^ ■\/3/2, 

lim Msn{z] A) / VA = kiis\n(kiix). (15) 



^ In Ref. [13, rather i>2 = is treated as an argument of 
Jacobi's elliptic funetion. In this notation Eq. mOI I becomes 
M(,y/v^sn{kz, U2). 
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FIG. 2. (color online), (a) The amplitude of the dynamical 
mass and an effective order parameter mavc in the solitonic 
phase as a function of 772 . (b) The magnitude of wave vector 
q = 2n/ip as a function of 772 in the solitonic phase. 



In the latter limit, va goes zero so that the magnitude 
of condensate |M(z)| vanishes while its form becomes 
sinusoidal. These two limits are given by exactly the 
same forms as we assumed in the previous section. We 
show in Fig. [2^a) the amplitude of mass function VAkA 
as a function of 772- Also depicted together is an effective 
order parameter 7Tiave defined by the root-mean-square 
mass averaged over the elliptic modulation period £p — 
AK{vA)/kA with K being the complete elliptic integral 
of the first kind: 



mavc = V{M{z)^) = kA ^1-E{va)/K{va). (16) 

E is the complete elliptic integral of the second kind. 
We also show in Fig. [2{b) the magnitude of wave vector 
q = 2iT/ip. From these figures, we clearly see that both 
phase transitions at 772 and 772^ are continuous ones. 

Finally we show in Fig. [3] the phase diagram in the 
[a2, sgn(a4)a|] plane which is essentially the same as the 
one displayed in figure 1 of Ref. [l^] ■ 
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FIG. 3. (color online). The phase diagram in the 
[q2, sgn(Q4)a|] plane. The horizontal axis represents 02 mea- 
sured in the unit [a^"^], namely 020^. For a4 > 0, the x^B 
phase (02 < 0) and the Wigner phase (02 > 0) are separated 
at 02 = by a second-order phase transition. For 04 < 0, 
there are three phases: the homogeneous chiral condensate 
indicated by xSB, the solitonic chiral phase marked by M{z), 
and the Wigner phase. Phase transitions separating these 
three phases, the dashed lines (black and red), are second 
order. For comparison, would-be first-order phase transition 
from the homogeneous xSB phase to the Wigner phase is also 
shown by the thin solid (gray) line inside the inhomogeneous 
phase. 



E. Singular behavior of thermodynamic quantities 
at the onset of soliton formation 



The phase transition at 772 = rf^ from the Wigner phase 
to the solitonic phase is second order, as clearly seen 
from the GL potential Eq. ([6]). What about the order 
of transition for a soliton formation at 772 = Tj\l Near 
the critical point v ^ \, the modulation period diverges 
logarithmically as Ip ~ 4/\(i/) ~ — 21og((l — i^)/8) so 
that the spatial separation between domain walls, a half 
period /2, also becomes large. The corresponding wave 
vector q — 27r/£p = ^kIv) behaves as 



5/6 

log((l-;/2)/16)- 



(17) 



We see that the wave vector drops to zero, inversely pro- 
portional to the logarithm of (1 — v), just as quark num- 
ber density in the GN model [2^ [29[ . The problem is 
how the elliptic modulus behaves as a function of 772 near 
the critical point. As 772 is approached, we see — > 1, 
fc^ 5/6 with A 4/5. From Eqs. and (fT2b)) . 

we can extract the expansion of v and fc^ in terms of 
(4/5 — A) and 772 — ?72- To first nontrivial order, we find 
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y'd - A) -f §(772-77!) 



(18) 



Also the relation between A and 772 — rj^ was numerically 
found to be linear in the vicinity of critical point (772 > 
77^), as A 9^ I - 2.406(7/2 - 77^). We obtain 



1 - 



73 '^iv^ 
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3V3 



111, 



(19) 



with a numerical value of kj = 2.38. Therefore, when the 
critical point 772 is approached from above, the inverse of 
the modulation length behaves as 



7rv/10/3 



log(§( 



V2 - m 



(20) 



We have checked to see that this analytical formula ex- 
cellently reproduces the numerical result displayed in 
Fig. m^b). This means that the separation of domain 
walls becomes logarithmically large while the wall thick- 
ness (the size of defect ~ 1/fc) does not change drastically 
near ri\. 

We now look at the singular behaviors in the thermo- 
dynamic quantities. These are determined by the behav- 
ior of derivatives of the GL free energy with respect to 
the parameter 7/2. The first derivative in the vicinity of 
the critical point can be computed up to the first few 
orders in the expansion in 772 as 



ws 



' A. _ 3^ 
12 4V/2 



_5_ 
12 



31og(^(r,2-»?i)) 



iv2 > vh) 



(21) 

We see that the first derivative is continuous. The second 
derivative is also worked out, and the first two dominant 
parts (either regular or singular) can be extracted as 



fiV2 



dv'i 



5+i 



3(r,2-r,i)log2 



i{V2-vl 



)) 



+ 



('72 < vl) 

{m > vl) 



(22) 

It is clear that the second derivative is not only discon- 
tinuous but diverges as the critical point is approached 
from the side of the inhomogeneous phase. The singular 
behavior is characterized mainly by the power law but 
with a logarithmic correction just as with the fermion 
number susceptibility in the GN model [2^ [2^ . We see 
that the transition is second order. We remark that the 
same singularity was obtained also in the NJL model [l^ 
where it was shown further that the vector-type interac- 
tion between quarks washes out the singular behavior. 

Before closing this section, we briefly discuss the origin 
of the singular behavior of the second derivative of ft. At 
T — the singularity is related to that of quark num- 
ber susceptibility [llllilii- To see this, we first note 
that at r = a net quark number (imbalance to qq con- 
densate) tends to accumulate in the wall surfaces where 
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the magnitude of the chiral condensate vanishes. Then 
considering the fact that the separation length between 
domain wahs diverges as £p/2 ^ 1/ log(7;2 — vl) while the 
size of the wall where quarks are present, 1/fc, stays al- 
most constant, the averaged quark number density may 
be approximately computed as 



(q^q) ~ no 



1/fc 



no 



£p + l/k log I 



(23) 



with no defined as the quark density in the absence of a 
condensate. Assuming uq does not change significantly 
in the vicinity of the critical point, the quark number 
susceptibility on the side of the inhomogcneous phase 
becomes 



d{q^q) di]2 



no 



d^i d/. (^2- r;!) [log (11(772-7,1))]' 



(24) 



This has exactly the same parametric dependence on r]2 
as Eq. ([m, suggesting that at T = the divergence of 
the second derivative of has originated in the divergent 
quark number susceptibility [l3|, i28|, i2^] . 



III. OTHER CRYSTALLIZATION PATTERNS 
IN ID AND HIGHER DIMENSIONAL 
MODULATIONS 

In this section we address the question of whether or 
not states other than the solitonic condensate, with ei- 
ther ID or higher dimensional modulations, are possi- 
ble. We work near the Lifshitz point, so we retain up 
to sixth order in GL expansion as in the previous sec- 
tion. In Sec. lnra we discuss the thermodynamics of the 
LO-like chiral condensate, Eq. ([5]), and the FF-like chiral 
spiral in some detail. Even though these two states 
are less favorable than the solitonic solution, the analysis 
still serves as an illustrative benchmark when we explore 
higher dimensional modulations. In Sec. IIII A[ we try a 
general Ansatz for the ID modulation, i.e., the conden- 
sate expanded in higher harmonics. For comparison, we 
also make a harmonic analysis on the solitonic state. In 
Sec. IIII Bl we introduce higher dimensional analogs of the 
LO-like chiral density wave, and see if such higher di- 
mensional chiral lattices can be realized near the Lifshitz 
point. 



A. Real chiral density wave and chiral spiral 

Here we discuss the thermodynamics of two typical chi- 
ral density waves: a real, LO-like sinusoidal chiral density 
wave of the form Eq. ([5]) and a FF-like chiral spiral char- 
acterized by a single plane wave 



(25) 



The imaginary part should be understood as the pseu- 
doscalar condensate, for one choice in the charge neutral 



channel: Mff(x) ~ —2G{{qq) + iiq^dir^q)). This is also 
called the "dual chiral density wave" abbreviated to the 
DCDW in the original paper In the following we re- 
fer to this state mainly as the FF state or chiral spiral. 
When the chiral condensate involves a finite imaginary 
part, we need to generalize the GL functional so as to 
allow it. The generalized functional for complex chiral 
condensates but with restriction to ID modulations was 
worked out first in the chiral GN model [13] ; the result- 
ing functional has terms asymmetric with respect to an 
interchange between the real and imaginary parts of the 
condensate, Re(77i) -f-)- lui(rn). These terms are respon- 
sible for the realization of the FF-type complex chiral 
condensate in the GN model ll^l and also in thie (1 + 1)- 

In 



m 3 



dimensional NJL model in the large N limit 
the NJL model in the l-f-3 dimension, however, these 
terms are shown to vanish p^ . Then the functional is 
cast in the following form in the chiral limit, using the 
same convention as Eq. (|4|): 



f |m|2 



, isgn(a4)(|m|4 + |mf ) 
6-K4|m|2|m'|2-HRe(m')2(m*)2 



(26) 

When the condensate is real [m(z)* = m(z)], this func- 
tional reduces to Eq. (jj]) with the restriction to the ID 
modulation, Af(x) — ^ M{z), substituted. As before, we 
concentrate on the case 0:4 < 0. In Fig.[4]wc show the am- 
plitude of masses (a), the magnitude of wave vectors (b), 
and corresponding thermodynamic potentials (c) for the 
FF and LO condensates, as a function of 772- The wave 
vector q is just k (ku) for the FF (LO) state, while for the 
solitonic state it is given by 27r/^p. From the figures, we 
see that the real LO phase is more favorable than the FF- 
like chiral spiral over the whole range of 7,2, while they are 
less favorable than the solitonic phase, which will often 
be denoted by the solitonic (SN) state hereafter. The FF 
and LO (and SN) states become degenerate at the onset 
of condensate at rj^ = 3/8. Near the point it is easy 
to perform the expansion in Afo. The result for the LO 
phase is already given by Eq. and that for the FF 
phase is obtained at the same order in AIq as [l3| 



n = (w(AfFF(x)))ws 



' 22 
\ 2 



A/2 



2^"0 



0(A/6). 



(27) 



The coefficient of A/q for the FF phase is twice as large as 
that of the LO phase, from which we see that the energy 
for the FF phase is higher than that for the LO phase at 
least near the onset. This may be reasonably understood 
along the same line as [s^ \^ where it is stressed that 
additional terms are required in the effective theory in 
order to recover the broken time reversal symmetry in 
the FF state, which results in an extra energy cost. 

Since the LO phase is more favorable than the FF 
phase, we concentrate on the LO phase in the follow- 
ing. The transition between the homogeneous x^B and 
LO phases is of first order. We can obtain the critical 772 
by solving miuM w(Arconst.) = minM,fc(w(Afo sinfcz))ws- 
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FIG. 4. (color online). (a) Effective order parameters, 
(b) wave vectors, and (c) potentials as a function of 772 for 
the LO (dashed, red), FF (dotted, blue) and SN (solid) states. 
The quantities corresponding to the homogeneous condensate 
are shown by dot-dashed lines. The solid lines are the same 
as those depicted in Fig. (2] 



This was worked out numerically and the location of the 
critical point was found as 



miLO O xSB) = 0.1519, 



(28) 



which is larger than r/l ^ 0.1389 for a soliton forma- 
tion, reflecting the fact that the LO phase is metastable 
against the SN phase. Crossing the critical point from 



the xSB phase to the LO phase, the magnitude of the 
condensate drops by about 20%. Also the ratio of the 
amplitude of mass to the magnitude of the wave vector 
in the LO phase just at the critical point has been found 
to be about 0.8. These can be summarized as 



Mo(LO) 



M,. 



t.(xSB) 



0.81, 



Afo(LO) 
fc(LO) 



0.80. 



(29) 



We note that all these ratios are the universal constants 
associated with the first-order phase transition between 
the xSB and LO phases, at the sixth order we are work- 
ing. 

We also remark that the first Ansatz of coimting Mq 
and k as the same order in magnitude is consistent at 
this transition. At the onset of the condensate, however, 
the amplitude Mg vanishes while k remains finite so that 
the derivative terms are more important than the bulk 
homogeneous terms. In contrast, at the onset of the do- 
main wall formation, the derivative terms play a minor 
role. 



B. Most general condensate with higher harmonics 

We here try to see if Jacobi's elliptic function is the 
most favorable solution among ID modulation patterns. 
We first note that Jacobi's elliptic function just gives 
a one parameter subgroup of solutions to the EL equa- 
tion. It should be stressed that the original EL equation 
dHl) is a nonlinear fourth-order differential equation while 
Eq. pT|) is basically the sum of three differential equa- 
tions each of which can be obtained from the second order 
differential equation which the elliptic function obeys. It 
is thus not obvious that it really covers all the solutions 
to the original EL equation. Keeping this in mind, we try 
the most general assumption for the spatially modulated 
chiral condensate. We set the condensate in the form of 
a harmonic expansion series as 

AfHH(^;{^H,MH})= M^sin{{2TT/tu)nz), 

n=l,3,5,--- 

(30) 

where {M„} and £h are variational parameters. We note 
that the even components are absent if we restrict the 
condensate to a half antiperiodic period M{z + ^h/2) = 
—M{z). This is quite reasonable for the case of the chiral 
limit. In fact we have checked numerically that the even 
components are vanishing even if included. Hereafter, we 
abbreviate the state characterized by Mhh as the higher 
harmonic (HH) state. We took into account up to the 
fifth harmonics in the expansion above and minimized the 
thermodynamic potential with respect to Af;P,M|^,M|^ 
and the modulation period ^h- 

The resulting energy of the HH state appeared to be 
very close to that of the SN state. In order to illustrate 
the fact, we depict in Fig.[5lja) the effective order param- 
eter mave = \/ (A/hh/sn(2;)^) as a function of 772, both 
for the SN and HH states: the solid line for the SN state. 
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and the dashed (magenta) hne for the HH state. We have 
labeled four representative points, 772 = 0.3749, 0.2569, 
0.1536, and 0.1389, by letters A, B, C, and D. In the 
figure, the coordinate (772, ^^^avo) at each point for the SN 
(HH) state is marked by the lower (upper) triangle. We 
see that the difference between the HH and SN states is 
visible only at point D which is very close to rji for the 
domain wall onset. 

The spatial profile of the HH-condensate, Mhh(z), was 
also observed to converge to that of the solitonic state, 
Msn{z); the difference is actually invisible already at this 
level of truncation in the biggest range of 772. It only 
becomes slightly manifest in the very vicinity of 772 ~ 
r}\ where the condensate becomes a collection of domain 
wall solitons. In Figs. OJb) and (c), we show the spatial 
profiles of the mass function at points A, C and D. Only 
at D do we show both Mhh and Msn, since at other 
points the differences are invisible. 

In order to see most efficiently how the HH and SN 
states diverge as the soliton onset is approached, we now 
make the harmonic analysis on the solitonic condensate 
Msn. First we note that the decomposition of the elliptic 
function in the Fourier-Lambert series is given as 



Ms„(z;{j/, k}) 



= 1,3,- 



M. 



i/, k) sin 
47rQ(z^)"/2 



27r77Z 



£p(i.,fc) i-g(z.)"' 



(31) 

where £p{i>,k) is the elliptic (real) period given by an 
explicit form £p(i/. A;) = AK{v)/k^ and Q{v) is the nome 
defined by Q{v)=exp{-iTK'{v)/K{v)) with K{y) and 
iK'{v) = iK{y/l — v"^) being the quarter periods of Jaco- 
bian elliptic functions. We compare {M^{v, k), £p{iy, k)} 
computed from Ms„(z) with {M^Ja}- In Table. HI we 
tabulate the values of potential, elliptic modulus, and 
harmonic coefficients computed from the SN state Msn 
and those for the HH state at four points introduced 
above. A, B, C, and D. From the table, we see that 
the energies of these two condensates are nearly degen- 
erate, actually the difference is visible only at the point 
D located in the vicinity of the onset of soliton forma- 
tion. Nevertheless it is true that the energy for the soli- 
tonic SN state is always smaller than that for the HH 
state. This means that even if we start from the general 
Ansatz Mhh (2), the optimized spatial profile of Mhh(-2) 
approaches that of Msn{z). This fact suggests that the 
elliptic function gives the absolute ground state. We rec- 
ognize that the closer to the point 7/2 we get, the more 
relevant the higher harmonic components become. Near 
7^2, Mhh (2^) starts to fail to approximate Msn (2); the de- 
parture of Mhh (2^) from A/sn(^) simply refiects the fact 
that truncated higher harmonic components (n > 5) can 
no longer be neglected. This is what actually is seen 
in Fig. [51^b) and (c); Mhh (2^) at points A and C weU 
approximates Msn (2), while it underestimates both the 
amplitude and the period at point D. 
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FIG. 5. (color online), (a) Spatial profile of the condensate 
Msn{z) at points A (772 = 0.3759; solid) and C (772 = 0.1536; 
dashed). For each case, the difference from A1uh{z) is invisible 
so it is suppressed from being depicted. The unit of horizontal 
(vertical) axis z (M) is set to the elliptic period (amplitude) 
of Msn(z), i.e., lp{r]2) {k{rj2)u{-rj2)) at each point, (b) Msn 
(solid) and Mhh (dashed) at point D (772 = 0.1389). The 
scales for M and the spatial coordinate z are normalized by 
the amplitude and the period of Msn{z) as in (a). 



C. Multidimensional modulation 

We now explore the possibility of realization of higher 
dimensional crystals favored over the ID solitonic state. 
We remark that the EL equation ([9]) was obtained after 
the restriction of condensate to one varying in ID. The 
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TABLE 1. Comparison between Ms-niz) and Mnni^z). The energy densities fisn and J^hh, the values of Fourier strengths 
{Af^'''^", M^''"", M^^''"}, and modulation periods 1-r and £p(i', fc) are tabulated at four representative points rj2 ~ 
0.375, 0.257, 0.154, and 0.1393 which are assigned by the letters A, B, C, and D, respectively. 





V 










Mr 


A/H 


Mr 


Mf 


Mr 


Mf 


A (r?2 = 0.3749) 


0.009 


5.13 


5.13 


-9 X 10""' 


-9 X 10""^ 


0.0111 


0.0111 


6 X 10"** 


6 X 10"* 


0(10"") 


0(10"") 


B {ri2 = 0.2569) 


0.458 


5.85 


5.85 


-0.0038 


-0.0038 


0.5291 


0.5291 


0.0077 


0.0077 


0.0001 


0.0001 


C {ri2 = 0.1536) 


0.821 


8.22 


8.21 


-0.0161 


-0.0161 


0.8673 


0.8665 


0.0563 


0.0559 


0.0039 


0.0038 


D (772 = 0.1389) 


0.991 


14.98 


10.48 


-0.0193 


-0.0192 


1.0672 


0.9781 


0.1951 


0.1077 


0.0455 


0.0132 



original EL equation in 3D is 
= A2a/(x) + 3AA/(x) - 10 [M{VMf + Af ^AM] 
+6r/2A/ - 6A/3 + 6^/^ 

(32) 

where the operator A + dy + is the three- 

dimensional Laplacian. This is a fourth-order nonlin- 
ear partial differential equation. The solution space for 
the partial differential equation is much wider than that 
for the ordinary differential equation. Therefore it is not 
trivial that Jacobi's elliptic function stays as the most fa- 
vorable structure when the restriction to ID modulations 
is taken away. Here we do not search for the formal solu- 
tion, but only try some specific crystals having the simple 
square or cubic symmetry. In order to demonstrate how 
the dimensionality of a crystal structure affects the free 
energy, we concentrate on the LO-type phase which has 
a simpler form than the elliptic function. We set multi- 
dimensional LO-type real condensates as 

Mlo;1d(x) = y2A'/osin(fcz), 

^^lo;2d(x) = A/o(sin(fca;) -I- sin(fcy)), 

-^^L0;3d(x) y^Afo(sin(fca;) -I- sin(fcy) + sin(fc2;)). 

(33) 

M^,0:2T) (x) is equivalent to the "egg-carton" Ansatz in 
|34| . Each condensate is characterized by two real pa- 
rameters Mq and q whose values are to be determined 
via minimization of 57. In Fig. EJa) and (b) we show 
effective order parameters and free energies. One can 
see that the energy is an increasing function of the di- 
mensionality. Also it is notable that the critical points 
at which the transitions from the Wigner to the crystal 
phases take place are common among all three states. To 
see this analytically, we expand in powers of Mq as 

^^LO;lD = (f -^)A.f2 + AM4 + 0(M6), 

r!LO;2D = (f -^)A'/2 + J^A/4 + 0(A./6), (34) 

f^LO;3D = (f -^)Af2+ l|A/4 + 0(A/6). 

Since the quadratic coefficients take exactly the same 
form, ^ — = ^(772 — ?72^), the critical point is shared. 
On the other hand, from the coefficients of quartic terms, 
we see that the free energies are on the order of the di- 
mensionality of modulation. The above formulas for en- 
ergy density can be easily generalized to the case with 
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FIG. 6. (color online), (a) The effective order parameters 
for the ID-LO (dashed, red), 2D-L0 (dotted, green), 3D-L0 
(dot-dot-dashed, magenta), and the solitonic state (solid) as 
a function of 772. For comparison, the quantity for the FF 
phase is also depicted by the thin solid line (blue), (b) The 
thermodynamic potentials as a function of 772. 

the LO condensate in an arbitrary dimension, d, defined 
by 

A/LO;dD(x) = y|A//o (sin(fca;i) + • • • + sin(fc.Td)) , (35) 

with X — (xi, X2, ■ ■ ■ , Xd) being the d-dimensional vector. 
In this case the thermodynamic potential looks like 

f^LO;dD = (f -^)A./2-fM^Mo^ 

, 5(4d^-16d+ll) a ^"^^^ 
48d5 ^"O- 



11 



From the dependence of the quartic coefficient on the di- 
mensionahty, (2d — 1)/Ad, we see that the energy is an 
increasing function of d near r]2 ~ rf^ ■ In fact, after min- 
imizing over Mq, the extremal of energy density becomes 



LO;dD — 



16) 



(37) 

Note that the higher dimensional LO phases are less fa- 
vorable than the ID one, yet they have smaller energies 
than the ID-FF (chiral spiral) state. 

We remark that the ordering of energies with the 
dimensionality is not universal, but rather specific to 
the LO-type. To see this fact, we first define the d- 
dimensional FF-type complex condensate as 



A'/, 



FF;dD 



(x) = 



(38) 



Then after working out similar calculations as in the LO 
case, we obtain the potential 



FF:dD 



^ 2 



3 ^ 7i.r2 I 1 71 r4 4d^ + 10rf-7 j r6 
Yg j JWq + 2^'-'0 1232 JWq ■ 



(39) 

In this case we see that the deference first appears at sixth 
order. Moreover, the coefficient (7 — lOd — M^)/{12d'^) 
is no longer a monotonic function of d; it has the ab- 
solute minimum at d = 1.4. Above this value, the en- 
ergy is a monotonically increasing function. Looking up 
the values (4^^ + lOd - 7)/(12d2) {0.55,0.60,0.58} for 
d ~ {1,2,3}, we see that the energy density is on the 

order of riFF;2D < f^FF;lD < ^^FF;3D- 

We also tried several other, more exotic Ansdtze for 
the higher dimensional crystal, including the 2D hexag- 
onal lattices both in the FF type and the LO type, but 
could not find the structure with energy less than the ID 
solitonic state. 



IV. PHASE STRUCTURE OFF THE 
TRICRITICAL POINT 

We now address the question of how the phase struc- 
ture changes going away from the tricritical point. As 
long as we restrict ourselves up to sixth order in the GL 
potential, the phase structure in the [a2, sgn(a4)Q;4] plane 
is simple such that any phase boundary is expressed by 
a line starting from the origin. There is no chance for 
these lines to meet again away from the origin. However 
this may not be the case in a realistic situation as stated 
in the Introduction. In this section we thus try to extend 
the previous analysis going beyond a minimal (sixth or- 
der) GL. and see if such nontrivial phase structure indeed 
shows up off the tricritical point. 

General eighth-order terms can be constructed by dif- 
ferentiating the sixth-order terms twice, and collect- 
ing terms which are unique up to total derivatives. 
At order ©(V^), only M''(VM)2 comes out, and at 
order ©(V) three terms [{\/ , M {\/ Mf I^M , and 
M2(AAf)2] show up. Finally at order 0{\/^) only the 



term (VAM)^ is possible. Thus the most general eighth- 
order terms to be appended to the GL Lagrangian Eq. ^ 
are collected as 



,5w8(M(x)) 



"86 



o 8 

+f^(VM)4 + —M{VMfAM 
8 8 

8 8 

(40) 

In NJL-typc models, only the quark loops contribute to 
the GL coefficients. The evaluation of them is tough but 
straightforward to work out. At the end of the day we 
arrive at the linear relations (up to total derivatives); 



(asf,, asc, asd, "Se, "s/) 



'U -1 18 14 l\ 
\^ 5 ' 5 ' 5 ' 5/ "8- 



(41) 



Thus it is only ag that we need to add as a new indepen- 
dent GL parameter. 

Here we briefly repeat the dimensional and scaling 
analyses in this case. Since ag ha-s dimension A~^ and 
is always positive, we use it for the energy scale. Then 
by use of jael we can again tune the magnitude of the 
sixth- and eighth-order coefficients. To be precise let us 
introduce the dimensionless GL couplings {772, '74} and 
the physical quantities {w, m, x} via 



a2 



ri"6i'i 


a 


1 

s . 


^. 




1/2 

"8 




1/2 



M = m 



(42) 



Then in these dimensionless bases, we have 



w = ^772^^ + |;774 ^m'* + (Vm)^ 

+ (^^6 ^ 5^2(y^)2 ^ 1 (A,„)2^ 

+ l{m^ + 14m4(Vm)2 - i(VTO)4 

+ ^m(VTO)2Am + ^m^iAm)^ + i(VAm)' 

(43) 

In this case, we have two independent dimensionless 
couplings (772, 774 )• When ag > the limit ag — >■ 0+ 
can be taken, and in this limit the thermodynamics 
should not depend on ag so that the physical quanti- 
ties depend on (772,7/4) only through the combination of 
772/774 = a2C(6/<^'i- This is the realization of (Q;2/a4) scal- 
ing in the absence of eighth-order terms that results in 
a linear relation between 012 and at any phase tran- 
sition. This translates into a linear relation between 772 
and ?74 in the limit of vanishing ag. We now expect a 
nonlinear relation between 77I and 772 due to the effect of 
ag by going away from the origin. From now on we throw 
away tildes on physical quantities, so we should keep in 
mind that they are dimensionless and properly rescaled 
quantities. 
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FIG. 7. (color online). Phase diagram for aa > 0: The phase 
boundaries in the absence of eighth-order terms (see Fig. [S]) 
are also repeated by thin dashed lines, for comparison. 



A. Phase diagram for ae, > 

We begin the discussion with the case of ag > 
which has a close relation with the previous analyses with 
ag = 0. The phase diagram for this case is displayed in 
Fig. [71 The phase boundaries are depicted by dashed 
lines. All the transitions are of second order. Just for 
comparison, we also show the phase boundaries obtained 
for the case ag = with two thin dashed lines, which cor- 
respond to two lines in Fig. [H As expected, we see that 
two critical lines are no longer linear in the plane, which is 
a simple consequence of the violation of the (772/^74) scal- 
ing. We can derive the analytic formulas for two critical 
lines. Analyzing the condition of the soliton formation 
min{7\/o,fe} /(Mo, fc, 772, m) 0''' with / defined by 



f{Mo,k,1]2,1l4) 



rf(^(M,n(z))) 
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FIG. 8. Phase diagrams for ag < 0: The horizontal (ver- 
tical) axis is Q2 (04) in the appropriate units shown in the 
parentheses, which is actually equivalent to the dimensionless 
coupling, ri2 {rji)- The phase diagram obtained with restric- 
tion to the homogeneous condensate (upper panel) and that 
computed without limitation to the homogeneous condensate 
(lower panel). Solid lines refer to first-order phase transitions, 
while dashed lines represent the second-order phase transi- 
tions. 



we easily arrive at the analytical expression for the soliton 
onset in a closed compact form, which can be expanded 
in powers of 774 < as 



V2 



36 



125 
1512 



3125 
28224 
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By use of original couplings, this translates into 
125^8 



0:2 ^6 



36 



al 



1512a2 
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3125ai 4 
28224^"" 



(45) 



(46) 



from which we clearly see that the eighth-order coeffi- 
cient as is responsible for the nonlinearity between a2 
and al . In the same way, we can also derive the compact 
formula for the critical line dividing the Wigner and in- 
homogencous phases. This can be expanded for 774 < 
as 



27 
160 



^74 



2187 
12800 



^74 



(47) 



Again the terms equal to or higher than the cubic term 
are caused by nonvanishing ag. We found that, despite 
these modifications of two critical lines, the qualitative 
feature of the phase diagram remains unchanged. In 
particular, the inhomogeneous phase is still dominated 
by the one-dimensional modulation characterized by the 
elliptic function Msn{z) = MQVsn[kz, v) with k ~ Mq. 



B. Going to new regime: Phase diagram for aa < 

We now have a closer look at the case ag < which 
is not connected in any limits of the previous analysis 
with Qfg = 0. In Fig. [51 we show the phase diagrams for 
ag < with and without restriction to the homogeneous 
condensate. 

The upper panel shows the phase diagram computed 
with restriction to the homogeneous condensate. In this 
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case, we have first-order chiral restoration even for a4 > 
because of negative ae- Moreover when 014 > 9/2 there 
are two types of x^B phases, one with a larger chiral con- 
densate, and the other with a smaller condensate which 
we label as the "XSB2" phase. The phase transition from 
the XSB2 phase to the Wigner phase is second order. On 
the other hand, the transition between the x^B and the 
XSB2 phases is first order. The magnitude of the chiral 
condensate jumps across the first-order phase transition. 
The gap between two condensates at the transition de- 
creases with increasing 014, and it eventually ends at the 
critical point marked by CP. The coordinate of this crit- 
ical point is found as (?72,?74) = 1/3). Accord- 
ingly, CP scales with ag as 
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(48) 



from which it is clear that point CP goes to the ori- 
gin (0,0) when ag — > 0^; it smoothly continues to the 
Lifshitz point in the region ag > 0. At point CP, the 
thermodynamic potential takes the following form in the 
dimensionless units: 



1 
648 



1 



(49) 



We see that the sixth-order derivative of di with respect 
to the order parameter m vanishes at this point. We 
expect large thermal fluctuation effects in the vicinity of 
CP. 

The lower panel of the figure shows the phase diagram 
where the limitation to the homogeneous condensate is 
taken away. A major part (but not all) of the first-order 
phase transition splits into two transition lines, and an 
inhomogencous phase shows up in between which is dom- 
inated again by one-dimensional modulation character- 
ized by the elliptic function. Point CP remains unaf- 
fected by the inclusion of an inhomogeneous condensate. 

In addition to CP, three points show up in the phase 
diagram. Let us closely investigate these points in the 
following. In the upper panel of Fig. [31 we present the 
magnified figure of the phase diagram. First, we notice 
that the Lifshitz point [L) moved to ^4 > from the ori- 
gin. The situation is the same at the point in which three 
phases, the homogeneous broken phase, the inhomoge- 
neous phase, and the Wigner phase, meet up. However, 
it now has the branch of the first-order transition sep- 
arating the inhomogeneous phase and the homogeneous 
XSB2 phase. The exact location of point L can be found 
analytically as (r/2,r/4) = (0,5/18) that means 



5 

L : (02,^4) = ( 0, — — 
' 18 as 



(50) 



Point (L) also will become degenerate with the origin as 
ag — > 0~; it extends to the Lifshitz point in the region 
ag > 0. 




FIG. 9. (color online). (Upper panel) A magnified version 
of the phase diagram in the (0:2, Q4) plane for ae < in the 
lower panel of Fig. [H] (Lower panel) The order parameter 
{M{z)^), and q — 2tt/£p as a function of 02 along the line 
04, — 0.273 which is displayed in the upper phase diagram by 
the dashed (magenta) line. 



There is also the critical end point (E) on the criti- 
cal line dividing the inhomogeneous phase and the ho- 
mogeneous broken phase. At this point the second- 
order phase transition turns into the first-order one. 
The location of this point can be numerically found as 
im,m) = (0.014,0.16), that is. 



E : (02,04) = 0.014 



0.16- 



(51) 



The most intriguing point that appears for ag < is 
the triple point denoted by T in the figure. At this point, 
three first-order phase transitions meet at once. The lo- 
cation is found numerically as (772, ^74) = (—0.016,0.27) 
which corresponds to 



T : (02,04) 



-0.016 



0.27^ 

as 



(52) 



Three different forms of the chiral phase compete with 
each other and coexist at this point. Any dimension- 
less ratio approaches to some constant when point T is 
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approached. We found numerically the following values 
M(xSB2) 



0.5 



lim 

(v2.v^}^T M(xSB) 



im^v^^T M(xSB) 
lim — 



0.369- •• , 

0.308- •• , 
5.001- •• , 



(53) 



where q = 27r/£p with i-p being the elliptic period. All 
these ratios are universal in the sense that they do not 
depend on any of the GL couplings. Thus we would say 
these ratios are associated with the triple point ( T) it- 
self. In other words, these values characterize the point 
(T). Therefore the values obtained here are model inde- 
pendent and should be shared by all theories which allow 
the existence of such triple point. 

All four points, CP, T, L, and E, discussed above get 
shrunk to the origin when the limit ag — ^ 0~ is taken, and 
they continue to the Lifshitz point in the region ag > 0. 

In the lowe r panel o f Fig. [5J we show the effective order 
parameter {M{z)-^), and g as a function of a2 along 
the section 774 = 0.273 which is depicted by a horizontal 
dashed (magenta) line in the phase diagram displayed in 
the upper panel. We see that the chiral restoration in this 
case proceeds via three steps: first from the xSB to the 
XSB2 phase, second from the XSB2 to the inhomogeneous 
phase, and last from the inhomogeneous phase to the 
Wigner phase through second-order phase transition. 



C. How does the GL map onto the (/i, T) phase 
diagram? 

The microscopic evaluation of GL couplin gs w ithin the 
NJL model can be summarized in the form [l0| 



(2^ 



(54) 



where cOm = 7rT(2m -I- 1) is the fermionic Matsubara fre- 
quency. The summation over the frequency can be done, 
resulting in the following expressions for a2 and 

a,^^- m^Nf J ^l^ME^I^-ME±tl, (55a) 



a4 



dp 1-/f(p-/j)-/f(p+m) 



(55b) 



dp fF{p-P-) + f'F{P+P-) 



(27r)' 



We can derive similar expressions for higher order co- 
efficients. The momentum integral in a2 {c^i) has 
quadratic (logarithmic) divergence so that we need a cut- 
off A. As a consequence, their functional forms become 
a-i = A2/i(^/A,T/A,.g) and = /2(/^/A,T/A) with 
g = Gh? as the dimensionless NJL coupling and {/«} 
being some dimensionless functions. Higher order co- 
efficients, a2n {n > 3), no longer have the UV diver- 
gence so that their functional forms can be summarized 
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FIG. 10. (color online). Mapping of the conditions of van- 
ishing GL coefficients to the (/i, r)-phase space. Solid (red) 
lines correspond to Q2 = for several difi'erent values of the 
NJL coupling; g = GA^, 2.0 (light), 2.5 (medium), and 3.5 
(heavy). The dot-dashed (blue) line presents the condition 
0:4 = 0. The intersection of the lines, Q2 = and 04 — 0, 
provides the candidate of the critical point. Area shadings 
are made according to the signs of Qe and ag. 
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FIG. 11. Conjectured, schematic phase diagram. The Lifshitz 
point splits into four critical points: the Lifshitz point, the 
critical point, the end point, and the triple point, labeled by 
L, CP, E, and T, respectively. 



as a2„(Ai,T) = /„(/i/T)/^2(„-2)^ (^^^ jg ^ function 
having three (two) arguments, while /„ for (n > 3) has 
only a single argument, /i/T. Thus the condition of van- 
ishing Q!2,i for n > 3 turns into a straight line in the (/i, T) 
plane if fn{x) = has a solution for a; > 0. 

It is easy to evaluate the vanishing conditions of 
0^2,4.6,87 ^iid the results are depicted in the (/i/A,T/A) 
plane in Fig. 1101 Solid (red) curves show the vanishing 
of a2 for g ~ 2.0 (light), g = 2.5 (medium), and g = 3.5 
(heavy), respectively The dot-dashed (blue) curve rep- 
resents the vanishing of 04. The intersection of these two 
curves, where a2 = 04 = 0, presents the candidate of the 
location of the critical point, in the absence of inhomo- 
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gcneous condensates. We see that the location shifts to 
higher temperatures with increasing the strength of the 
NJL couphng, g, while its /i coordinate shows nonmono- 
tonic dependence on g. 

The long-dashed (green) line and the dashed (magenta) 
line show the lines for vanishing ag and respectively. 
Accordingly, the condition (as > 0, ag < 0) is realized in 
the shaded area, while ag changes its sign in the lightly 
shaded area where (ag < 0, < 0) is satisfied. The 
region without shading corresponds to (ag > 0, ag < 
0). The present analysis suggests that if the "would-be" 
critical point is located in the shaded area, this point 
may split to the four independent points, CP, T, L, and 
E discussed in the previous section. If this is the case, 
the corresponding (//, T)-phase diagram might become 
the one schematically depicted in Fig. [TT] 

We remark that this interesting situation is realized for 
very large values of the NJL coupling. The NJL coupling 
empirically has taken ranges from g = 1.8 to g = 2.5 de- 
pending on the presence/absence of the 't Hooft coupling 
0, M, SH 113 ■ Then it would be unlikely to have the triple 
point in the standard version of the NJL model. It would 
be interesting to investigate if the situation could be re- 
alized in the extended version of the NJL models such 
as one incorporated with the Polyakov-loop [43l - l47| . the 
other type interactions, the diquark degrees of freedom, 
etc. 



V. SUMMARY 

We have made an extensive analysis on inhomogeneous 
chiral condensates using the Ginzburg-Landau approach. 
We first reviewed the work by Nickel [Ifll with special em- 
phasis put on the mechanism of the formation of a chiral 
domain wall. We have also discussed singular behaviors 
of the thermodynamic quantities in the vicinity of soliton 
formation. 

We then analyzed the most general condensate of the 
ID modulation, Mjm^z) expanded in harmonics. We 
confirmed that AfuH seems to converge to the solitonic 
state characterized by Jacobi's elliptic function by in- 
creasing the number of harmonics. It well approximates 
the solitonic condensate except for in the vicinity of a do- 
main wall onset, already at the level of truncation with 
the first few harmonics. This strongly suggests that the 
elliptical solitonic state is the most favorable structure 
within ID modulations. 

We then explored the possibility of realization of mul- 
tidimensional chiral crystals. Taking several Ansdtze for 
multidimensional structures, we compared free energies. 
We have shown that, in the case of real condensates of 
the LO-typc, the free energy density is an increasing func- 
tion of the dimensionality of modulation d. On the other 
hand, it has been confirmed that the energy takes a min- 
imum for d = 2 in the case of complex condensate of 
the FF-type. We derived analytic expressions for the 
effective potentials, both for the real and complex con- 



densates forming a simple square and cubic lattices. The 
order of energy was correctly understood based on the 
GL expansion near the critical point. 

Since we have restricted our analyses to simple lattice 
structures, the quest for such crystal structures remains 
the subject of continuing interest. In particular, because 
simple ID structures are known to be unstable against 
the thermal fluctuation [l^ , we should clarify how the ID 
structure could be stabilized if we really have no suitable 
multidimensional structure. 

We have investigated the phase structure away from 
the tricritical point by extending the previous GL func- 
tional to one expanded on the order parameter and its 
spatial derivative up to eighth order. With this exten- 
sion, we analyzed phase structures both for > and 
ae < 0. For ag > 0, we have seen that two critical lines 
enclosing the inhomogeneous phase change their slopes 
in the [a2, sgn(Q;4)a4] plane by going away from the Lif- 
shitz point. This is a natural consequence of including 
the eighth-order terms. We have derived the analytical 
expressions for the two critical curves. However no new 
phase structure was found to appear even largely depart- 
ing from the Lifshitz point. 

On the other hand, once ag changes its sign to neg- 
ative, ag < 0, several qualitative changes appear im- 
mediately in the (012, Q!4)-phase diagram. We have seen 
that the Lifshitz point located at the origin for ag > 
splits into four critical points; these are the critical point 
[CP), the triple point (T), the end point (E), and the 
new Lifshitz point (L), which has a branch of the first- 
order phase transition. The distances among these points 
grow with increasing the magnitude of |Q;g|. There are 
two kinds of homogeneously symmetry-broken phases in 
this case: one with a larger chiral condensate xSB, and 
the other with a smaller chiral condensate which we la- 
bel by xSB2. They are separated by a first-order phase 
transition extending from the triple point (T) to the 
critical point (CP) at which the phase transition ter- 
minates. Most intriguing is the emergence of the triple 
point (T) where three different forms of the chiral phase, 
xSB, XSB2 and the inhomogeneous (solitonic) phase, are 
present. Any dimensionless ratios of physical quantities 
approach some specific values when the triple point (T) is 
approached. These values are universal, model indepen- 
dent, being associated with the triple point itself. We 
have extracted numerically the universal ratios associ- 
ated with the triple point (T). 

There are several possible extensions of the current 
work. We have restricted our analyses to the chiral limit. 
It is straightforward to extend our work so as to include 
the effect of a finite current quark mass along the line of 
[30| . It may also be interesting to try the mo st g eneral 
complex condensate, the twisted chiral crystal |3l| . even 
though a simple chiral spiral condensate of the FF type 
turned out to be the most favorable structure within a 
ID-NJL model [s^l- Also the interplay between the chiral 
condensate and color superconductivity is of particular 
interest since it was shown that there is an entanglement 
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between them via the axial anomaly in QCD [49[. Its 
effects have been proven to be so drastic that it brings 
a rich phase structure near the critical point (49l - l52| . A 
more interesting, but challenging issue is the first prin- 
ciple determination of the Ginzburg-Landau coefficients 
with the lattice QCD. Extension of the present work in 
these directions deserves future work. 



J. Wambach arc acknowledged for their interests in this 
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members of the theoretical quark/hadron group at Tokyo 
University of Science for useful conversations. 
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